Dark solitons in a two-component Bose— Einstein condensate 
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The creation and interaction of dark solitons in a two-component Bose-Einstein condensate is 
investigated. For a miscible case, the interaction of dark solitons in different components is studied. 
Various possible scenarios are presented, including the formation of a soliton-soliton bound pair. 
We also analyze the soliton propagation in the presence of domains, and show that a dark soliton 
can be transferred from one component to the other at the domain wall when it exceeds a critical 
velocity. For lower velocities multiple reflections within the domain are observed, where the soliton 
is evaporated and accelerated after each reflection until it finally escapes from the domain. 



The realization of Bose-Einstein Condensation (BEC) 
in weakly interacting atomic gases [0 has opened the 
possibility to investigate nonlinear properties of atomic 
matter waves. In this respect, several remarkable results 
have been reported such as the experimental observation 
of four-wave mixing in BEC 0], and the realization of 
vortices |,[§, and dark solitons in BEC §,§. 

A dark soliton in BEC is a macroscopic excitation of 
the condensate with a corresponding positive scattering 
length, which is characterized by a local density minimum 
and a sharp phase gradient of the wavefunction at the po- 
sition of the minimum . The shape of the dip does not 
change due to the balance between kinetic energy and 
repulsive atom-atom collisions. The recent experimental 
creation of dark solitons in BEC by means of the phase 
imprinting technique (^||, has posed several fundamen- 
tal questions concerning the dynamics, stability and dis- 
sipation in such systems [p|-[lO||. Also, the interaction of 
two solitons in a BEC has been experimentally addressed 

In the recent years, the development of trapping tech- 
niques has allowed the creation of multi-component con- 
densates. These are formed by trapping atoms in dif- 
ferent internal (electronic) states JlJ,[ll|. The multicom- 
ponent BEC, far from being a trivial extension of the 
single-component one, presents novel and fundamentally 
different scenarios for its ground-state wavefunction 
and excitations |l5f| . In particular, it has been experi- 
mentally observed that the BEC can reach an equilib- 
rium state characterized by the phase separation of the 
species in different domains ]l3| . 

In the present Letter, we analyze the creation, prop- 
agation and interaction of dark solitons in a two- 
component condensate using analytical and numerical 
methods. In this more complex scenario, novel phenom- 
ena can be expected, as it has been already reported in 
the context of nonlinear optics H . We show that the dy- 
namics of the soliton interaction is completely different 
from the single-component case. In particular, two dark 
solitons in a single-component BEC always repel each 
other fig], whereas the opposite is true for two solitons 



interacting in the two-component case |17| . We consider 
two situations. For the first one, the two components are 
miscible, and one soliton is created in each component. 
We show, both analytically and numerically, that this 
system presents nontrivial dynamics, which includes the 
formation of a soliton-soliton bound state. The proper- 
ties of such a binary bound state, in particular the pe- 
riod of its oscillations, should be sensible to dissipation 
effects and therefore constitutes an excellent tool to an- 
alyze these effects. In addition, the dissipation can be 
studied in a much clearer and controllable way than in 
current experiments with single-component BEC, since 
the soliton-soliton interaction keeps the solitons within 
a central region of the trap of several healing lengths. 
In the second scenario, the propagation of a soliton in 
a two-component BEC with domains is considered. In 
particular, we demonstrate that the soliton can be trans- 
ferred from one component to the other at the domain 
wall. Below a critical velocity the soliton is reflected, per- 
forming multiple oscillations inside of the corresponding 
domain. At each reflection the soliton is partially evap- 
orated by emission of phonons, and accelerates until it 
eventually escapes though the wall. 

In the following we consider a trapped BEC with two 
components, where the dynamics takes place only in one 
dimension due to the strong trap confinement in the 
transverse direction. This approximation is valid if the 
mean-field interaction is smaller than the typical energy 
separation E± in the other directions jl8| , and has been 
successfully employed in the analysis of dark solitons in 
single-component condensates M. For sufficiently low 
temperatures the dynamics is well described by two cou- 
pled Gross-Pitaevskii equations 
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where V(x) is the trap potential, /ij is the chemical po- 
tential of component j, gji — AirTi 2 aji/mS is the cou- 
pling constant between the components j and I with the 
transversal area S, m is the atomic mass and dji the 
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scattering length between j and I (j, I = 1,2). 

We consider the situation in which a dark soliton is 
created in one or in both components. To generate dark 
solitons we use the well established method of phase im- 
printing [EJj| , although other methods |Q could in prin- 
ciple be employed. This method consists of applying a 
homogeneous potential U generated by the dipole po- 
tential of a far detuned laser beam to one part of the 
condensate wavefunction of one component. The poten- 
tial is pulsed on for a time t p , such that the wavefunction 
locally acquires an additional phase factor exp(iUt p /h). 
The pulse duration is chosen to be short compared to 
the minimal correlation time of the system h/fj,, with 
fx = max(/ii, fi2)- This ensures that the effect of the light 
pulse is mainly a change in the phase of the conden- 
sate, whereas changes of the density during this time can 
be neglected. Note, however, that due to the imprinted 
phase at larger times an adjustment of the density in the 
condensate appears, leading to the formation of the dark 
soliton and also additional structures. In principle the 
described method can be employed to selectively create 
a soliton in only one of the components, or to create two 
different solitons in each one. The creation of a soliton in 
one component modifies the density in the other one, due 
to the coupling in Eq. (EJ) . But as we show below, under 
appropriate conditions, solitons can also be created in a 
two-component condensate. 

We consider first the case in which one soliton is cre- 
ated in each component. This situation can be analyti- 
cally studied by employing a variational approach |^l|,^2| . 
We assume for simplicity the solitons are moving in a ho- 
mogeneous condensate of densities ni = 71% — uq. This 
situation corresponds to a very elongated trap with equal 
concentrations of both components. Additionally we as- 
sume that the coupling constants are gu = g\2 = gi2 — 
g. The latter assumption matches well the experimen- 
tal conditions |||n|. As variational wavefunctions the 
single-component soliton solutions are considered 
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where fi\ = fX2 = A*, g = dq/dt, 2q(t) denotes the relative 
distance between the solitons, c s = y/ 'gno/m the sound 
velocity and Iq = %/ ^/gnom is the coherence length for a 
single component. Eqs. ( pa| ) and (^) represents a kink- 
antikink situation, i.e. when the phase fronts of the soli- 
tons are facing each other. The case kink-kink, when 
both phase fronts are in the same direction, is also dis- 
cussed below. The previous expressions describe a sym- 
metric situation around x — 0, or equivalently they de- 
scribe the system in the center of mass frame. 

The problem of solving Eq. (Eh can be restated as a 



variational problem [|21|,[22| , corresponding to the station- 
ary point of the action related to the Lagrangian density 
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Our goal is to find the equation which governs the evolu- 
tion of q(t). In order to do that, we insert the variational 
Ansatz (2a),(pb|) into Eq. (||), and calculate an effective 
Lagrangian L = J dxC which becomes 
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where b(q,q) = \f\ — q 2 1 'c 2 a {q / 'Iq) . The equation for q(t) 
is provided by the Euler-Lagrange equation ^ (§§) ~~ 
= 0, which we have solved for different initial condi- 
tions, obtaining the trajectories in the phase space (q, q) 
(Fig. EJ). For small deviations from (0, 0), the system be- 
haves periodically , i.e. the solitons form a bound pair 
(soliton molecule). On the other hand the free trajecto- 
ries are characterized by the acceleration of the approach- 
ing solitons, and the deceleration of the outgoing ones. 
For q > 2lo, L becomes the sum of two single-soliton 
Lagrangians and q vanishes. For practical purposes the 
trajectories can be considered periodic if they cross q = 
at q < 21q, and free otherwise. The free trajectories close 
to the periodic ones, are squeezed together at q = and 
q = q c = 0.73c s , which constitutes the critical escape 
velocity. 

The kink- kink situation is however different. In this 
case the solitons are propagating in the same direction 
with velocities q\ and q2, such that in our formalism 
2? = <?2 — Si- For this case L, which has already been de- 
rived in the analysis of the stability of optical vector dark 
solitons plfl , takes the form of the first line of Eq. (|J). 
For small velocities also bound soliton solutions appear, 
oscillating for small deviations from (q = 0, q = 0) with 
a frequency ^/8/15/z/TL However, contrary to the kink- 
antikink case, the solitons can never break the molecule. 
This is reflected in our formalism by the appearance of 
a singularity at the escape velocity q = c s /^/2, which 
is the half of the sound velocity for the homogeneous 
two-component gas c Si 2 = V2c s . The singularity reflects 
the fact that the solitons cannot move in the labora- 
tory frame faster than c Sj 2- The special kink- kink case 
(q = 0, q = 0) coincides with the optical vector dark 
soliton solution [E_7 23 . 

We have studied the creation of a soliton in each com- 
ponent using phase imprinting where different initial con- 
ditions corresponding to different phase imprintings are 
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obtained. We consider the case of a BEC with equal 
number of atoms in both components, TV = 10 5 , in a 
box trap and with gu = 922 = 1.05<7i2. This choice 
allows a homogeneous region of equal densities for both 
components, and therefore provides a better quantitative 
comparison with the analytical model. However, for more 
general non-homogeneous situations a similar qualitative 
picture has been observed in simulations. 

We restrict ourselves here to the case kink-antikink. 
The creation of a soliton in one component perturbs the 
density of the other one. We observe that if the phase 
imprinting is applied in order to create two solitons which 
are initially separated by distances larger than approxi- 
mately 4?o, the fluctuations in the densities prevent the 
formation of the solitons. Therefore, the solitons have to 
be created initially with q(Q) < 21q. The soliton in each 
component induces a local density increasement centered 
at the position of the soliton in the other component. 
In other words, the solitons are filled by the other com- 
ponent. Therefore, they become wider and slower than 
a corresponding single-component soliton. This fact in- 
troduces some quantitative corrections to the analytical 
estimates, although such corrections are in fact small. 

Fig. H shows the case of an initial q — 0.8 lc s . For such 
velocity, the solitons move periodically around q = 0, 
i.e. they are forming a soliton molecule, as described 
above. The separation of the solitons, which can reach 
4/im, depends on the initial velocity given by the phase 
imprinting. We have numerically found a critical velocity 
q c = 0.83c s at which the solitons become free. This veloc- 
ity is in good agreement with the variational approach. 

Fig. U shows the evolution of the solitons correspond- 
ing to an initial velocity q = 0.89c s (dashed line in Fig. 
|3|). The solitons indeed move apart much slower with 
q = 0.095c s . The closer q is to q c the larger the soli- 
ton deceleration is in agreement with our analytical re- 
sults. Therefore, the deceleration is indeed an effect of 
the soliton-soliton interaction and not a consequence of 
the filling of the soliton by the other component. We have 
also depicted In Fig. || the trajectory after the reflection 
from the box boundaries, in order to illustrate the be- 
havior when both solitons collide. As predicted from our 
variational approach, it can be observed that the solitons 
are accelerated when approaching each other, and decel- 
erated after crossing. The maximal velocity at q = is 
comparable to the critical velocity. 

In the last part of this Letter we analyze numerically 
and analytically the soliton propagation in separate do- 
mains. If the relation between the coupling constants 
and densities is appropriately chosen, separate domains 
of each component can be created We consider the 
case in which a soliton is created in one of the components 
and move towards the domain wall. In order to illustrate 
the different possible scenarios, we study the situation in 
which <7i2/.9ii = 1.7, giijgw = 0.96 for different initially 
imprinted velocities. Both components have equal num- 



ber of atoms N = 10 5 . A sufficiently fast soliton will 
be transferred through the domain wall into the other 
component. However, if the velocity is sufficiently low, 
the soliton is reflected at the domain wall, as shown in 
Fig. § This fi gure shows the case of a box trap with 
the initial soliton velocity q — 0.15c s . At each reflection 
the soliton is partially evaporated in the form of phonons 
in the second component. The latter induces an accel- 
eration until the soliton eventually escapes the domain. 
The critical escape velocity can be estimated from sim- 
ple energetic considerations, assuming that the soliton 
must overcome a potential barrier induced by the sec- 
ond component at the domain wall. This gives a critical 
velocity q t = \J (gn — gn)/4Slo, where S is the transver- 
sal area. In the considered example, the analytical value 
qt = 0.19c s is in excellent agreement with the numerical 
one qt = 0.16c s . When the soliton is transferred a back 
action of the soliton on the domain is observed. This 
introduces density fluctuations and perturbations in the 
domain walls, which slightly modify the critical velocity. 
The latter can produce a retrapping of the soliton in the 
original domain, as observed in Fig. |J. 

In this Letter we have shown the rich behavior of soli- 
tons in two-component BEC. The two components pro- 
vide solutions such as bound solitons and the possibility 
to create extremely slowly moving ones. We have ana- 
lytically studied the dynamics of the system with a vari- 
ational approach, and determined the possible scenarios. 
We have finally analyzed a two-component BEC which 
contains domains, and showed that depending on the 
physical parameters a dark soliton can be either trans- 
ferred or reflected at the domain wall. 

Several interesting problems remain, however, open. 
Among them, we stress especially two. In the present 
Letter we have analyzed a ID system. If the ID con- 
ditions are not strictly fulfilled, dynamical instability is 
expected IJhJ. In the new scenario with two-component 
condensates the properties of such instability should be 
altered. A second interesting problem is given by the 
dissipation of the oscillatory motion. The two solitons 
radiate phonons when oscillating. Contrary to the case 
of other binary systems, the radiation will increase the 
elongation of the oscillations, until eventually breaking 
the soliton molecule, and therefore these systems could 
be an excellent probe for the dissipation effects [M . 

We should finally stress that the effects considered 
here appear for realistic situations and can be experimen- 
tally analyzed with the state of the art technology. The 
creation of dark solitons constitutes a well established 
technique for the case of a single-component BEC. We 
have numerically simulated the phase imprinting mech- 
anism in a two-component BEC and demonstrated that 
this technique can also be applied in that situation [ B5j . 
Since the solitons are indeed wider due to the presence of 
the second component, some of the predicted effects, as 
for example the appearance of a critical escape velocity, 
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could be experimentally observed in a non-destructive 
way. Others, however, as for example the soliton oscil- 
lations, could require the opening of the trap, and sub- 
sequent condensate expansion. The dynamics of such 
expansion will be the subject of a separate investigation. 
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gemeinschaft (SFB 407), TMR ERBXTCT-96-002, and 
ESF PESC BEC2000+. Discussions with J. Anglin, J. 
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H. Schmaljohann and K. Sengstock are acknowledged. 



due to experimental limitations, there are alternative 
ways of creating soliton molecules. E.g. the same phase 
imprinting in both components, together with a displace- 
ment of one of the components, also leads to the same 
effect. 



FIG. 1. Phase map of the kink-antikink relative motion. 
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FIG. 2. Density of component 1 for the kink-antikink case 
with g(0) = and q(0) = 0.81c a . Darker regions are those 
with less density. Component 2 is the mirror image of com- 
ponent 1 around x = 0. 



FIG. 3. (Left) Density of component 1 for the 
kink-antikink case with q(0) = and q(0) = 0.89c s . Darker 
regions are those with less density. Component 2 is the mirror 
image of component 1 around x = 0. The dashed line is the 
soliton trajectory in a single-component condensate. (Right) 
Detail of the collision region. 



FIG. 4. Interaction with a domain wall of a soliton initially 
created in component 1 with q — 0.15c,,. 
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